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Abstract 

The main issue in this paper is to discuss a parity property that 
appears in the expressions for the distorted spectrum of the thermal 
Sunyaev-Zel'dovich effect. When using the convolution integrals method 
involving scattering laws we argue that the distorted spectrum contains a 
new term, which is an odd power of the frequency. Such a term, absent in 
the conventional approaches, implies a crossover frequency which differs 
in value form the ones reported in the literature by a significant, in prin- 
ciple observable, amount. Also, such term casts doubt on the demanding 
need of computing complicated relativistic calculations. The relationship 
of our approach with the existing calculations is discussed. 



1 Introduction 

The Sunyaev-Zel'dovich (SZ) effect is a signature left in the cosmic microwave back- 
ground (CMB) by collapsed structures The thermal Sunyaev-Zeldovich effect 
arises from the frequency shift of CMB photons that are scattered by the hot electrons 
contained in structures such as galaxy clusters. The frequency dependence of this 
effect results, for a given line of sight, in an intensity decrease in the Rayleigh- Jeans 
region of the CMB spectrum and to an intensity increase at Wien's region. The ef- 
fect has been detected in observations of some rich, X-ray luminous clusters. Besides 
this overall effect in the light spectrum, of most interest is the value of the crossover 
frequency, which is an observable quantity commonly used in the determination of 
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cosmological parameters [3] HJ. This observable is a strong ingredient in our foregoing 
discussion. 

The study of the non-relativistic and relativistic thermal effect using analytical 
methods has led many authors to agree, one way or another, upon the fact that the 
distorted spectra arising form Compton scattering can be expressed as a power series 
containing terms which are even powers of the frequency v, except for a first order in v 
that is present due to the diffusion approximation El El- The parameter z = 
used as a discriminant to indicate when the relativistic corrections are important is 
usually kept up to second order effects. In this paper, we challenge these results. We 
claim that if the correct basic physics behind the inverse Compton scattering is used, 
a term which is odd in v, in fact of order z 2 v 3 must appear, leading to a couple of 
singular results. The first one is that the crossover frequency, v c , is subject to a 2.8 
percent correction with respect to the value predicted by diffusion, in comparison to the 
1.19 percent correction obtained with the relativistic equation. These values are taken 
for fcT e = 5KeV. This could be presumably tested with precision observations. The 
second implication is that, even when kT e — 15KeV , this odd power term, not only 
favors a value of v c with an 8.4 correction with respect to the 4.3 percent predicted by 
relativistic corrections, but clearly indicates that the curves for the distorted spectra 
are practically indistinguishable from those obtained by numerical methods, specially 
in Wien's region. Moreover, it is precisely in this limit where claims have been made 
pointing out the opposite In short, we sustain the old claim issued by Sunyaev 
25 years ago |1U|. stating that there is no need to perform relativistic corrections to 
obtain the correct spectra even at the upmost values for kT e = 15KeV. 

These arguments will be subsequently developed in sections two and three of this 
paper. 



2 Gaussian scattering laws 

We begin this paper by recalling the approach followed previously to discuss the SZ 
effect in terms of a scattering law |1U| |11| . This method is based upon the fact that 
the distorted radiation spectrum which originates from inverse Compton scattering 
between electrons and photons in a certain cluster is given by: 

I(u)= / I {v)G 3 {D,v)dv (1) 



Here, I{v) is the scattered radiation off the plasma, I {v) = 2h ji (exp{yp-) — 1) 1 the 
undistorted spectrum, where Tr is the CMB temperature, and G s (v, v) the scattering 
law. For the thermal non-relativistic effect we have argued that the form Gs{v,v) is 
given by JO] 

G a {p,v) = {l-T)5{p-v)+TG[p,v) (2) 
where r is the optical depth for a thin plasma, and 

1 f P-(l-2,)„ .2 

G(p,u)= e v^<„) . ( 3 ) 

\/2-ko(v) 

2 lz T 2 2 

Here, a{v) = = 2zV is the square of the width of the spectral line at 

frequency v, and m e and T e are the mass and temperature of an electron and the 
electron gas, respectively. 
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Figure 1: CMB distortion for the case of Eq. (4) (dashed line) compared with the 
relativistic curve, Eq. (5) (solid line), here kT e — 7.5KeV . Both curves are practically 
identical deep in the Wien region. The frequency is given in Hz and — is given in 



The first term in Eq. Q is self explanatory, it represents the probability that a 
photon traverses the plasma unscattered. The second term is a modified Gaussian 
probability function describing the scattering of a photon with an incoming frequency 
v by an electron whose average kinetic energy is \kT e . The reason of why the outgo- 
ing frequency is shifted by a factor 2zv arises precisely form the fact that the average 
frequency of a photon scattered by electrons with speeds u = /3c will exhibit a tem- 
perature blue shift given by 2zv |10| 

Inserting Eqs. J2J and © into |T). and upon integration one is lead to the 
result that 



AIM dl„ 2 d 2 I„ o22 d 2 I 2 3 d 3 I 2 2 z/d 4 / 
r av ov 2 av 2 ov s 2 w 4 

to second order in z and first order in r. Eq. |0J is a useful expression in order to 
discuss the SZ thermal effect. Nevertheless, several characteristics of this result must 
be underlined. Here, of course, AI(u) = I(v) — I {u). 

Setting aside the fact that to first order in z Eq. JIJ is identical to the diffusion 
approximation described by the Kompaneets equation m |5] |12|. a fact that has been 
analyzed in Ref. |13|. we have other interesting features. The first one is the presence 
of a third order derivative term in the right hand side of Eq. and the other one 
is that the last term leads to the somewhat surprising fact that the Wien side of 
the distortion curve is practically identical to its relativistic counterpart This is 
clearly exhibited in Fig.(l) for kT e = 7.5KeV. 
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A careful examination of the analytical expression for the relativistic thermal SZ 
effect derived by Sazonov and Sunyaev shows the absence of a third order derivative 
term. Indeed, Eq. (12) of Ref. [Sj for a static cluster (V = 0) is identical to the 
expression: 



^ = (-2 z + f z*) v «g. + ( z - %z>) , 2 ^ 

7 J„4 8 4 7 D W 



a fact that can be easily proved by cumbersome but straightforward algebra. Before 
attempting to analyze if a third order derivative term should arise in the distortion 
curve, let us go back one step and ask how is it possible to obtain an equation for 
AJ ^ without the third order derivative term. We don't question the first derivative 
term, since it appears already in the diffusion approximation. The non-existence of 
the i/ 3 term can be achieved by writing the scattering law for G s as follows: 

G 3 = (1 - t)5{v - (1 - 2zt)v) + tG'{v, v) (6) 

where 

G'{p,v) = — ^ e " ( 7^ )2 (7) 
V2T0- 



Physically, this implies that G' '(&, v) is a purely Gaussian function peaked at v — v 
|14| and that the central limit theorem governs the scattering of photons by electrons 
|15| . Nevertheless, the equation would also imply a blue shift in the frequency for the 
unscattered photons. In fact, any even function G (a) in the variable a = would 
suppress the odd order higher derivative terms, as in Eq. ©. Although we believe 
that this is physically incorrect, substitution of Eqs. © and 10 into Eq. Q leads to 
an interesting result, namely, 

AJ(«/) a dl 2 d 2 I 1 a 4 9 4 7 

= ~ 2zl/ ^T + zv 77 +o zv "a-T ( 8 ) 

r w Oi^ 2 2 

In Fig. (2) we exhibit the effects of this result compared with the Kompaneets 
approximation and with the full relativistic curve of Sazonov and Sunyaev , Eq. 
for kT = WKeV. In Fig. (3) the curve arising form Eq. |gj is included. The reader 
may wonder why all this fuss about the higher order in z corrections but we believe 
that the underlying physics is basic for understanding the effect and, as we shall see 
below, also to grasp the nature of the relativistic corrections to G s ■ 



Thus, we are brought back to the original question, namely, if there is a third 
order derivative present in Eq. Q. To obtain a guide for replying we went back 
to a detailed examination of two papers, one by Stebbins written six years ago [S| 
and the publication of Sazonov and Sunyaev on relativistic corrections to both the 
thermal and kinematic SZ effects Although based on different approaches, both 
led to results which turn out to be identical for the thermal SZ effect. The same 
parity properties were found in related works [Sj - [Hj. Indeed, Eq.(10) of Ref. is 
essentially an expression for the distorted spectra written in terms of the occupation 

number n 7 = (e kT R — 1) instead of the intensity noticing that 

^ = ^ (9) 
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Figure 2: CMB distortion for the case of Eq. (8) (short dashed line) compared with 
the Kompaneets approximation (long dashed line) and the relativistic curve, Eq. (5) 
(solid line), kT e = WKeV. The order z 2 curves are practically identical deep in Wien's 
region. The frequency is given in Hz and ^ is given in units of 2( ^«) . 
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Figure 3: The same as is Fig. (2) now including a comparison with Eq. (4), (dot- 
ted line). The third order derivative term induces an increase in the crossover 
frequency and a better agreement with the relativistic curve at high frequencies. 
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where r is the optical depth, identifying e = hv and assuming that hv <C mc 2 in 
that equation. The expansion of all terms of Eq.(lO) of Ref. 6 is a laborious but 
straightforward calculation that shows that the third order derivative in I disappears 
and that the resulting expression is identical to Eq.JSJ. The surprising result of the 
comparison between the non-relativistic limit of Eq.@ and Eq. J1J is that the third 
order derivative in I ought to be absent. In our context, this would imply that, for 
low temperature clusters kT e ~ 5KeV, Eq. (SJ should be essentially correct, which 
we object because it is physically questionable. 

This leads us to a different test, namely the value of the critical frequency v c defined 
by AI(v c ) = 0. One gets, by the Kompaneets equation that v c — 220.027 GHz, a value 
which according to the relativistic correction increases, for kT e = 5KeV, only to 222.67 
GHz, a 1.19 percent correction. However, if the third order derivative is present in 
Eq.JU, as we assert, v c = 226.39 GHz, giving a more significant correction of 2.81 
percent. The simple physics behind Eq.lQJ is quite attractive but merely theoretical 
up to this point. Notice also that from the observation of the results, it is really hard 
to see which is in agreement with observations, all follow the same pattern, specially 
in Wien's region. 

Finally, a comment on the last term of Eq.JlJ. Our method, which is non- 
relativistic, yields a numerical factor of instead of the jg factor obtained in refer- 
ences and using explicit relativistic corrections. The question is if this difference 
and possible improvements can be obtained using the scattering law approach in the 
relativistic case. We turn to this question in the following section. 



3 Relativistic scattering law 

From the discussion of the previous section it is rather clear that in order to introduce 
relativistic corrections into the scattering Kernel G a (v, v) two facts must be considered. 
First, the energy of an electron has to be written in its relativistic form E — mc 2 ^— 1) 
and second, their velocity distribution is no longer given by Maxwell's distribution, 
but by its relativistic extension. Nevertheless, since the relativistic parameter z is 
small even at high temperature clusters, fcT e = 15.1KeV, we may simply generalize 
the scattering law, Eqs. © - ©, allowing for z 2 corrections for the peak shift and 
a(v), namely 

G r( v, „) = (!- r)8ip - „) + -^Hi^^L ) (10) 

In this case, the evaluation of the expression I(v) = I {v)G r {v,v)dv leads to a 
generic Kompaneets type equation. This is achieved performing the change of variable: 

y/2(r(v) 

with normalization factor 



and identifying 



N= H{a)da (12) 

J — oo 

Au(a) = V2a(v)a - ei> (13) 
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Figure 4: CMB distortion for the case of Eq. (4) (solid line) compared with Eq. 
(5) (short dashed line) and with the Kompaneets approximation (long dashed line) 
, kT e — Vb.lKeV . The crossover frequency would have a large increase due to the 
presence of the third order derivative term. The frequency is given in Hz and — is 

given in units of 2 . 

The resulting distortion curve reads: 

= ___ ev §± + ,28\ + (^)i^ f°° a 2 H(a)da 
+ ( — ) sI^J-oo a H{a)da+ { m )-gj£ J_ oc a H{a)da 

Eq. is a particular case of Eq. 1141 for e = — 2z. a 2 (u) — 2zu 2 making H(u. u) — 
_J___ e n/2<t(^j Qf coursej the series can be trivially taken to arbitrary order 
in z n . It is interesting to notice that the coefficient of the third order derivative 
term in Ea ll4|l is independent of relativistic effects up to second order in z and would 
yield —2zv 3 even taking into account relativistic corrections. In Fig. (4) we exhibit 
the effects of Eq. gj compared with the Kompaneets approximation and with the 
full relativistic curve of Sazonov and Sunyaev for kT — 15.1KeV. The critical 
frequency shift correction will now increase from 4.347 to 8.414 percent, basically due 
to the presence of the third derivative term. 



4 Discussion 

The present paper not only shows that the v 3 term appearing in Eq. @ is essentially 
enough to account for the relativistic features of the SZ thermal effect without entering 
into long and sophisticated calculations, but provides a bridge between the convolution 
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integral approaches to the thermal SZ effect and the expansions in terms of derivatives 
of the undistorted number occupation . The calculation of convolution integrals in 
the logarithmic space, which has been used by other authors |3] prevents the 
establishment of expressions such as Eq. 11411 that may relate the approaches. On 
the other hand, the calculation of the convolution integrals in the frequency space fills 
this gap. Similar techniques applied to CMB distortions have been discussed in other 
contexts |16| . 

To end this work we wish to underline the fact that the presence of the third 
derivative term is directly directed to the frequency shift included in the kernel of Eq. 

If we take into account that, on the average, a scattering produces a slight mean 
change of photon energy [3] |1U| . then direct use of the central limit theorem leads to 
Eq. ©. 

This work has been supported by CONACyT (Mexico), project 41081 — F. 
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